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ABSTRACT 

One of the difficulties in designing controllers for infinite- dimensional systems arises from 
attempting to calculate a state for the system. In this paper it is shown that Galerkin 
type approximations can be used to design controllers which will perform as designed when 
implemented on the original infinite-dimensional system. No assumptions, other than those 
typically employed in numerical analysis, are made on the approximating scheme. 
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1. Introduction 


There are computational difficulties, apart from the theoretical problems, to designing 
controllers for a class of systems whose dynamics are described by partial differential equa- 
tions. Unlike finite- dimensioned linear dynamical systems, which are often easily solved for 
the state of the system, a closed form solution for the state of an infinite- dimensional system 
can be computed only in the simplest of situations. This is due to the difficulty of solving a 
partial differential equation. In general, it is necessary to uBe a numerical approximation to 
the state of the system when designing a controller. 

This leads immediately to several questions regarding robustness. Is the use of finite- 
dimensional approximations a valid technique when designing controllers for 
infinite-dimensional systems? If so, which approximation methods can be used to design 
controllers which will perform as designed when implemented on the actual system, and how 
high an order is required? 

A wide class of projection methods, which includes mode truncation and finite element 
Galerkin techniques, is considered in this paper. It is shown that, if the order of an approxi- 
mation to a system is sufficiently high, then the approximating system is stabilizable, and a 
controller designed using the approximation will stabilize the original system. Furthermore, 
the closed loop response of the system will be close to the the closed loop response of the 
approximating system. This result is independent of the technique used for controller design. 

Other researchers have studied convergence of solutions to the Riccati equation [eg., 9, 
10, 12], but have assumed either uniform stability or stabilizability of the approximating 
systems. Balas [1, 2] assumes convergence of a stabilizing feedback control. However, the 
assumptions on an approximation scheme required for the results in this paper to hold are 
not stronger than those typically imposed on any numerical technique used in computer 
analysis. The key assumption is that the original infinite-dimensional system is stabilizable. 

We consider semigroup control systems of the following form: 

x(t) = Ax(t) + Bu(t), 

y(t) = Cx(t) t (1.0) 

x(0) = x 0 ®o € D(A) C X. 

A satisfies the Hille-Yosida theorem on the Hilbert space X and hence generates a strongly 
continuous (C 0 )-semigroup T(t) on X. Control and observation are assumed to be bounded 

I.C., 

B e B(i2 ra ,X) and C 6 B(X,R>). 
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This formulation includes a wide class of systems; in particular many problems concerning 
vibrations in structures with point sensing and actuation can be placed within this frame- 
work. 

The inner product on X is indicated by ( , ) and the corresponding norm by, || ||. Inner 
products on different spaces are not distinguished as the context indicates which space is 
meant. We use the following definition of internal stability. 

Definition 1.1: The Co-semigroup T(t ) is a-stable if there exist constants M and a > a 
such that ||T(<)|| < Me~ at for all t > 0. 

The possibility of a > 0 is included to allow for situations where the purpose of con- 
troller design is for purposes other than simple stability (t.e., a = 0). For instance, while 
vibrating structures are generally open-loop stable, the settling time of these systems is often 
unacceptably large. 

Several definitions will be required. 

Definition 1.2: The pair (A, B ) is cr-stabilizable if there exists a bounded linear operator 
K : X —* R m such that A — BK generates a a-stable semigroup. 

Definition 1.3: The pair ( A,C ) is <j- detectable if there exists a bounded linear operator 
F : R v — *■ X such that A — FC generates a a- stable semigroup. 

Definition 1.4: The system (A,B,C) is jointly a- stabilizable/ detectable if ( A,B ) is 
a-stabilizable and (A, C) is a-detectable. 

Define the linear space A to consist of functions of the form 

OO 

h = h a (t) + hiS(t - U). (1.1) 

i=l 

where h a E Li[0, oo) and the sequence {hi} is absolutely summable. The norm of a function 
in A is 

ii*i- rn.wi*+f:M<°°- (i.2) 

J0 i= 0 

A SISO system is said to be Loo-stable or bounded-input bounded-output (BIBO) stable 
if and only if its impulse response belongs to the set A. Moreover, the gain of an BIBO-stable 
system is equal to the norm of its impulse response in the set A [6]. The symbol A is used 
to denote the set of Laplace transforms of distributions in A. 

Define M( A) to be the set of matrices with elements in A. A MIMO system is Loo-stable 
if and only if its impulse response belongs to M( A). The L^-gain of a m x p matrix H is 
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the induced norm defined by 


\\H\\a ^max^p^llA. 

’ 3=1 

Consider the common situation where the system is already stable and the aim of con- 
troller design is to improve its transient response. Let us specify a real number <7 > 0, which 
is the minimum acceptable stability margin. We define a system to be a — Loo-stable if and 
only if its impulse response f(t) is of the form (1.1), and instead of (1.2) satisfies the more 
stringent conditions 

/o° 

2^ |/ii|exp(at;) < oo, / \h a (t)ex.p(crt)\dt < oo. (1.3) 

i=0 JO 

Let A a denote the set of all distributions satisfying, and let A„ denote the set of Laplace 
transforms of distributions in A^. 

We now define precisely what is meant by external stabilization of a system. Suppose 
P is a given system, for which we wish to design a controller C , arranged in the familiar 
feedback configuration shown in Figure 1. It is readily seen that the 2x2 transfer matrix 
H(P, C ) which maps the pair (u 1; u 2 ) into the pair (ex, e 2 ) is given by 


H{P,C) = 


(/ + PC)- 1 -P(I + CP)- 1 ' 
C{I + PC)- 1 (/ + CP)- 1 


The feedback system, or alternatively the pair (P, C), is said to be stable if each of the four 
elements in the above matrix belongs to the set S of stable transfer functions. Definition 
of S depends upon the application. Thus the closed loop system is BIBO-cr-stable if and 
only if all four elements belong to A ff . It is important to note that the present definition of 
stability is symmetric in P and C. Thus P stabilizes C if and only if C stabilizes P. 


u2 



Figure 1. Feedback System 
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The problem is to design a controller for the system (1.0) which achieves such objectives as 
external stability. Since this generally requires knowledge of the system state, the difficulty 
in solving partial differential equation systems and many other semigroup systems means 
that an approximation to the state must be used. In the next section we define a class of 
approximation techniques known variously as Galerkin and projection techniques. 

2. Approximation Technique 

Suppose we have a sequence of finite-dimensional subspaces Xh with C X^ for 
hi > h 2 . We assume that Xh C D(A). When finite elements are used as a basis for the 
approximating subspace Xh, then A should be understood in its weak or extended sense [7, 
20]. Define P h x as the orthogonal projection of x E X onto the finite-dimensional subspace 
Xh- The parameter h converges to zero as the order of the approximation is increased. Since 
all norms on Xh are equivalent, we will assume that the norm on Xh is that inherited from 
X. Since Ph is an orthogonal projection, it follows that 

llft*ll < Ml- 

The approximating system is (Ah, Bh,Ch) where 

Bh PhB 

C h := <7| Zfc 

and Ah is the bounded linear operator defined by 

A h x h = PhAxh for all x h E X h - 

In common with [3], the operators Ah, Bh,Ch and the semigroup Tk(f) generated by Ah 
are only defined on Xh- 

We make the following assumptions: 

Al) The projection operators Ph converge strongly to the identity on the Hilbert space 
X. That is, for all x £ X 

lim II PhX — s|| = 0. 

h--*Q 

A2) The approximation scheme approximates the solution space not the operator A, in 
that x E Xh implies Ax E Xh so that 

AhX = Ax for x E Xh- 

Modal truncation and many applications of the finite element method satisfy this assumption. 
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A3) A core C of a closed operator A is a linear space contained in the domain of A with 
the property that the set of elements (x, Ax'), x £ C is dense in the graph G( A) of the 
operator A ([15] pg. 166). In other words, if we define an operator L to be the restriction of 
A to C so that 

D(L) = C 
Lx — Ax, x € C 

then A is the closure of L, or equivalently, G(A) = G(A). 

We will assume that there exists a core C for A such that 

lim ||P/,A® — AhPkX || = 0, for all x G C. (2.1) 

Such an approximation scheme is said to be consistent. 

We will further assume that the semigroups Th(t) generated by Ah are uniformly bounded, 
that is, there exist real numbers H,M > 1 and u> such that 

||T fc (t)|| < Me ut for all h < H. (2.2) 

(Uniform boundedness of the approximate semigroups is generally referred to as “stability” 
in the numerical analysis literature.) Consistency and uniform boundedness are required for 
convergence of the approximation [e.g., 16, 18, 20] i.e., for all e > 0, t > 0 and for all x 6 X 
there exists H such that 

\\PhT(t)x — T>i(t)P/ix|| < e for all r G [0,t] and h < H. 

Assumptions (A1)-(A3) are satisfied by typical approximation methods such as finite 
element techniques and modal truncation. 

In subsequent sections the original semigroup control system will be assumed stabilizable 
and/or detectable. However, no assumptions about the stabilizability or detectability of 
the approximations are made. In the next section it is shown that stabilizability of the 
approximating systems can be deduced from stabilizability of the original semigroup control 
system. 

3. Stabilizability of Approximating Systems 

We require the following theorem which guarantees the existence of a solution to the 
infinite- dimensional Riccati equation if the pair (A, B ) is stabilizable. 

Theorem 3.1: [17] Let A and B be as defined above and let L,Q be self-adjoint and 
positive definite bounded linear operators. If there exists a solution u* to the optimal control 
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problem with the quadratic cost function 

J(u):= [ [(Lx(t),x(t)) + {Qu(t),u(t))]dt (3.1) 

Jo 

which satisfies the condition that J{u*) be finite, then there exists an unique self-adjoint and 
non-negative solution II E L(X, X ) to the infinite- dimensional Riccati equation 

A*U + TlA-UBQ- l B*Tl + L = ® (3.2) 

where 0 is the zero operator. Let K = Q~ l B* II; then A — BK generates an exponentially 

stable C 0 - semigroup. D 

Clearly an admissible solution exists if and only if (A, B ) is stabilizable. 

Theorem 3.2: Assume we have a semigroup control system (2.0) and an approximation 
scheme satisfying assumptions (A1)-(A3) in Section 2. Assume that the system is stabiliz- 
able. Then, for sufficiently small h, the approximating systems are stabilizable t'.e., there 
exists an operator K h such that A h - B h K h generates an exponentially stable semigroup on 
X h . 

Proof: For each approximating subspace, define a self-adjoint operator 11^ : Xh —* Xh 
as 

Ilh := PhUPh 

where II is the solution to the Riccati equation (3.2). It will be shown that, for small enough 
h, 11^ solves the finite-dimensional Riccati equation 

A^II/i + II/iA/, — UhBhQ -1 BhUh + Lh = ®h 

where Q h is the zero operator in B (X h ,Xh) and L h is a self-adjoint positive definite operator 
(as yet undetermined) in B(AT/,, X^). 

Define the operator R(Lh ) : Xh — * Xh by 

R(L h ) = + n h Ah - n h B h Q- l B' h Tl h + L h 

where Lh £ B(A\, Xh) is as yet unspecified. Let us calculate a representation for R(Lh) 
with respect to some basis for Xh- Using (3.2), we have 

(*, R{L h )<t>i) = ( <f>i , Q(h)h) + (h, Lh<f>j) - Lh) 

where 

Q{h) := TiBQ-'B'Tl - P h HP h BQ- l B*P h UP h . 
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Now define Lh • Xh — > Xh as 


L h x h = Lx h - Q{h)x h for all x h £ X h . 

Then R(L h ) = 0 on B(X h ,X h ) since 

(&, R(Lh)<t>j) = 0 for all &, <f>j. 

In order to show that II/, solves a Riccati equation for (Ah, Bh ) for small enough h, 
it remains only to prove that Lh is self-adjoint and positive definite for small h. Now, 
Hm h >0 || (/ - P h )B || = 0 since B operates on a finite-dimensional space. It follows that [22, 
Theorem 4.14] 

lim||<2(h)|| = 0. 

The fact that Lh is self adjoint follows easily from the definition of Q(h). Since L is positive 
definite and self-adjoint we can choose H so that for all h < H 

(x h , LhXh) > C 2 (x h , X h ) 

for some c ^ 0. Thus Lh is self-adjoint and positive definite for all h ^ H , and II/, . Ph H Ph 
solves the Riccati equation for (Ah, Bh) with weights Q and Lh if h is sufficiently small. 
Defining Kh := <5 _1 ^II/,, Ah - BhKh generates a stable semigroup on Xh for h < H. □ 

Banks and Kunisch [3] derived a similar result for parabolic systems using different as- 
sumptions. Rosen [19] considered the problem of convergence of solutions to approximating 
Riccati equations. However he considers a much smaller class of problems: his proof applies 
to cases where the open loop system is stable and n is a Hilbert-Schmidt operator. Although 
the result given here is weaker in that convergence of H/, to n in an operator norm is not 
shown, the theorem is considerably more general. 

The following corollary is a consequence of the two results above. 

Corollary 3.2a: Assume that (A,B) is cr-stabilizable and let n be the solution to 
the Riccati equation for (A + cI,B) with weights Q and R as in equation (3.2). Define 
U h := PJlPh and let K h = Q~ l B* h Tl h . Then for sufficiently small h, A h - B h K h generates a 
a-exponentially stable semigroup on Xh- D 

The results above for cr-stabilizable systems have obvious extensions to c-detectable 
systems since (A, C) is c-detectable if and only if (A*,C*) is cr-stabilizable. 

Theorem 3.3: Let (A, B, C) be a semigroup control system and suppose we have an ap- 
proximation scheme which satisfies assumptions (A1)-(A3) of Section 2. Assume furthermore 
that (A,B) is cr-stabilizable. Choose n £ B(AT,X) such that n solves the Riccati equation 
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(3.2). Define K := Q~ l B*Tl and let 5(f) to be the semigroup generated by A — BK. Simi- 
larly, define Kh := Q~ x B^PhUPh and let 5/,(f) be the semigroup generated by Ah — BhKh- 
Then 5j,(t) converges to 5(f) uniformly on bounded intervals; i.e., for all e > 0,f > 0 and 
for all x 6 X there exists H such that 

||P/,5(t)x — Sh(r)PhX || < e for all r e [0,f] and h < H. 


Proof: 

Ilft-M < IIbq-'b'II ||n|| ; =c.- ~ 

Using (2.2) and [18, Theorem 3.1.1] it follows that there exist real numbers H, M > 1 and 
u> such that 

||5fc(f)|| < Me^ +Mc ^ for all t > 0, for all h < H . (3.3) 

Define Ao := A — BK and Aho := Ah — BhKh . Let C be the core for A such that the 
consistency condition (2.1) is satisfied. Since BK is a bounded operator, C is also a core for 
A a . The resolvent of (A I — A a ) is bounded for Re(A) > — a, and so (A/ — A 0 )C is dense in 
X for ReA > — a. It follows that for every x € C, we have 

lim || AhoPhX - P h A oX || < lim WPhBQ^B^PhUPh - II)x|| = 0. 

In order to show convergence of {5/,} to 5 it necessary to introduce a semigroup which is 
defined on all of X. Define Aho — AhoPh- Since Aho is a bounded linear operator it generates 
the C 0 -semigroup §h{ f) defined by 

= E p-it = E = ShWp,. 

k=0 K - k = 0 K - 

Since ||5j,(f)|| < ||5fc(<)||, 5fc(t) also satisfies (3.3). Now, for x E C, we have 

lim - i4x|| < lirn \\AhoPhX - PhA 0 x\\ + || P h A a x - A 0 x\\. 

The first term was shown above to tend to zero . The second term converges to zero since 
Ph converges strongly to the identity. It follows that 

lim sup ||5h(t)x — 5(t)x|| =0 for all x 6 X, for all r > 0 

h~*° 0 <t<r 

by the Trotter-Kato Approximation Theorem [18, Theorem 4.5]. Therefore, for all x € 
X,t >0, 

lim sup ||5^(t)PhS — Pfc5(t)x|| = 0. 

h-*0 o <t<r 

□ 
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The following theorem is now immediate. 

Theorem 3.4: Suppose that the system ( A , B, C ) is stabilizable and that assumptions 
(A1)-(A3) hold. Assume also that the approximation scheme also satisfies the additional 
condition (A4): 

Urn sup || T' h {t)P h x - P h T'{t)x || = 0. 

0<t<T 

Then for sufficiently small h, the Riccati equation 

Afcirfc + 7 VhAh — KhBhQ~ l B^Kh + PhLPh = 0/i (3-4) 

has a unique non-negative self-adjoint solution T*h • The sequence of operators tt^ converge 
strongly to the solution to the exact equation (3.2). Define the optimal feedback operators 
Kh ■= R^BfJIh and K := R~ 1 BTl) the semigroup generated by Ah — BhKh converges to 
the semigroup generated by A — BK . 

Proof: By assumption, hypotheses (H1)-(H2) in [12] are satisfied. Since Theorem 3.3 
implies uniform stabilizability/detectability it follows that the sequence of solutions to the 
approximate Riccati equations are bounded by some real M x [12]. Since {717,} 6 B(H,H ) 
the sequence contains a weakly convergent subsequence. It is clear that the limit of this 
subsequence is also self-adjoint and non-negative. In fact, the limit is the solution to the 
exact Riccati equation, II [Theorem 6.7, 8]. A proof identical to that of [Theorem 6.9, 8] 
yields the remainder of the result. □ 

The fourth assumption can easily be shown to hold for a wide class of problems. How- 
ever, while assumptions (A1)-(A3) are almost always satisfied by a numerical scheme which 
is satisfactory for approximation purposes, (A4) may not be satisfied. This can lead to 
unsatisfactory controller design. An example of this situation is discussed in [4]. 

4. Convergence in the Graph Topology 

The graph and gap topologies [22, 23] arose from a need to define convergence of possibly 
unbounded operators. The basic idea was outlined in Kato [15]. If M , N are closed linear 
subspaces of a Banach space X then £(M , N ) is the smallest number 6 such that 

8 > dist(u, N ) = inf ||u — u|| for all u 6 M. (4-1) 

The gap between M and N is defined by 

&(M,N) := ma x(£(M, N),6(N,M)). 
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If S and T are closed operators from X to Y then the gap between their graphs as closed 
subspaces of X x Y is well defined. The difficulty in defining a gap for unstable systems is 
that the plant may not be represented by a mapping between two Banach spaces. However, 
some bounded inputs u, will be mapped to bounded outputs y. Let U be the set of stable 
inputs and outputs and let P be the impulse response of a scalar plant . We will assume 
that U is a Banach space. The symbol denotes convolution. Define 

Du := {u 6 U such that y €E U} 


and the operator 

P\j * u := P * u, 

with domain Du- Then the graph of P\j is 

G(P\j ) = {(u, y), y = P * u, u € Du}. 

It is easy to show that G(P\j) is closed if P is a closed operator. The graph of a plant is 
now defined to be the set of bounded (in the appropriate sense) input-output pairs, i.e., 

G v (P) = G(Pv) = {(u,y)eV m +’, y = P * u}. 

In the case of BIBO stability, 

Define the gap between operators as 

d(P 1 ,P 2 ) = ^(Gu(Pi),G u (P 2 )). 

With this metric the space of closed operators from U to U, C(U, U) becomes a metric 
space. Convergence of a sequence P n — ► P is defined by d(P n ,P ) — > 0. If a sequence of 
plants {Pi} are stable, then convergence in norm to P is equivalent to convergence in gap. 
The topology induced by the metric d(Pj, P 2 ) is thus an appropriate generalization of the 
uniform topology to unstable plants. 

Note that the above definition of “stable graphs” includes the graphs of non-causal plants. 
A more useful topology for our purposes, which considers only linear causal plants, can be 
obtained through the frequency domain approach of Vidyasagar [22]. Let the set of stable 

A A A 

plants S be A, P a plant with a right coprime factorization ( N,D ) over A. Define G 00 (P) 
to be the Laplace transform of elements of Goo(P). Then 

G „{P) 


10 



or 

G s (P) = {(Dz,Nz),zeLZ}. 

The graph topology for unstable plants is defined in [22] using right coprime factorizations 
(r.c.f.’s). Essentially, a basic neighborhood in the graph topology of a plant with a r.c.f. 
(N 0 , D 0 ) is the set of all plants with a r.c.f. (N, D ) where N is close to N a in norm and 
D is close to D 0 ■ Define S to be the usual algebra of stable rational functions. For finite- 
dimensional plants the graph topology is the topology induced on the quotient space S/S 
by the set of stable functions, S. Further details can be found in [22]. Zhu [24] presents a 
generalization of the graph and gap topologies to subsets of general quotient fields. Here, 
we consider only a particular set. 

Callier and Desoer [5] define a quotient algebra of transfer functions of distributed systems 
which is useful for studying control of semigroup systems. 

We say that a transfer function / in A is bounded away from zero at infinity if for every 

sequence {s/t} with Re(s) > 0 and 

1**1 00 

we have 

lim inf |h(sjfe)[ > 0. 

k-*oo 

Define the sets 

A_ — {he A a for some a > 0} 

A_ i00 = {h G A_, / bounded away from zero at infinity} 
and the corresponding quotient algebra 

B = [/i — a/b, a G A„, b G A__ |0 o}- 

The assumption that b G A_ >00 implies that functions in B possess only a finite number of 
unstable poles. The extension to the multivariable case is straightforward, and we define 
M(B) to be the set of transfer function matrices with entries in B. 

An extension which restricts the region of stability to Re(s) < — a is also straightforward. 
The algebra A is replaced by A„ in the above definitions and 

B 7 — {h — a/b, a G A^—, b G Ar,oo}‘ 

Jacobson [13] and Jacobson and Nett [14] have shown that semigroup control systems 
which are <j-stabilizable and/or cr-detectable have transfer functions in M(B a ). This frame- 
work is thus appropriate for studying semigroup control systems using frequency domain 
methods. Furthermore, every function P in M( B„) has both a left- and a right-coprime 


11 



factorization over M(A.&). Thus Vidyasagar’s graph topology is well-defined for semigroup 
control systems providing they are <7-stabilizable and/or -detectable. 

Theorem 4.1: Suppose Fh is a sequence in M(B„), and that F € M(B a ). Then the 
following statements are equivalent. 

(i) { F h } converges to F in the graph topology. 

(ii) There exist a r.c.f ( N,D ) of F, and a sequence of r.c.f.’s (Nh,Dh) of F, \ such that 
N h — ► N and Dh — * D in 

(iii) There exist a l.c.f. ( N , D ) of F, and a sequence of l.c.f.’s (Nh, Dh) of Fh such that 
Nh—>N and Dh — + D in M( X a ) . 

Proof: Since all transfer functions in M(B a ) have left and right coprime factorizations 
[5], Vidyasagar’s proof [21] for rational functions that convergence in the graph topology is 
equivalent to convergence of coprime factorizations extends in a straightforward manner to 
transfer functions in M(B a ). □ 

The following result can be loosely paraphrased as follows: A family of plants Fh can be 
robustly stabilized by a compensator C which stabilizes some nominal plant F if and only 
if F h converges to F in the graph topology. Furthermore, the closed loop response of the 
feedback pair ( F h ,C ) converges to that of (F, C). In other words, the graph topology is the 
weakest topology in which feedback stability is robitst. 

Theorem 4.2: Let Fh be a sequence of plants in M(B„) . 

(i) Suppose F h converges to F e M(B„.) in the graph topology. Let C 6 M(B„) stabilize 
F. Then there exists an H such that C stabilizes Fh, for all h < H, and moreover, the 
closed loop transfer matrix H(Fh, C) converges to H(F,C ) in M(A. a ). 

(ii) Conversely, suppose there exists a C £ M(B a ) which stabilizes F h for all h < H, and 
that H(Fh, C) converges to H(F , C). Then converges to F in the graph topology. 

Proof: Identical to that in [21]. □ 

This has obvious implications for controller design using approximations. Failure of a 
sequence of approximations to converge in the graph topology implies that at least one of 
the following must occur (1) simultaneous stabilization of Fh and F is not possible for all h 
sufficiently small or (2) convergence of the closed loop response H(Fh, C ) does not occur. 
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Theorem 4.3: Let ( A,B,C ) be a <7-stabilizable/detectable semigroup control system, 
and suppose we have a sequence of approximations (Ah, Bh, Ch) satisfying assumptions (Al)- 
(A3) of Section 2. If the system is cr-stabilizable and C is any controller which stabilizes P, 
then for sufficiently small H, (P h , C ) is cr-stable for all h < H. Furthermore, the closed loop 
response of (Ph, C) converges in norm to that of (P, C ). 

Proof: The theorem will be proven by using Theorem 3.3 to show that the approximate 
systems (A*, Bh, Ch) converge to (A, B, C ) in the graph topology, and then applying Theorem 
4.2. 

Let feedback operators K and Kh be as defined above so that Aho := Ah — BhKh generates 
a a-stable semigroup Sh(t ) for sufficiently small h and Ao := A — BK generates a a-stable 
semigroup S(t). For small h, right <7-coprime factorizations for the approximate systems are 
given by 

N h (s) = C h R(s : Aho)B h , 

D h (s) = I - K h R(s : Aho)B h 

where R(s : A^) indicates the resolvent of Aho [21]. A right coprime factorization for the 
original system is [13] 

N(s ) = CR(s : A 0 )B, 

D(s) = 1- KR(s : Ao)B. 

The result follows if it can be shown that Nh(t ) — ► N(t) and Dh(i) —* D(t ) in the norm 
on A a . Convergence of the numerators is proven first. 

Since {Sfc(t)} converge uniformly on bounded intervals to S(t), there exist H and M > 1 


such that for some 7 > a 



||S(t)|| < Me' 7 * for all t > 0, 


and 

11^(011 ^ Afe -7 * for all t > 0, for all h < H. 

(4.2) 

Therefore, for 

any e > 0 there exists t such that 



jf exp(<ri)||Af(() - W k (()||di < 

(4.3) 

Also, 

sup ||W(t) - JV k (i)|| < sup ||C|| ||S(t)B - S k (()B||. 

0<*<T 0<t<T 

(4.4) 


It follows from Theorem 3.3 that this also approaches zero. Combining statements (4.3) and 
(4.4) we see that 

feoij ex P( crt )H iV ( t ) “ ^(011* = 0 
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and so Nh converges to iV on A». 

The proof showing convergence of the denominators is somewhat lengthier and requires 
use of a lemma in [16]. Defining 

F h (t) := ( P h U-U)S h (t)P h B , 


we have 

D h (t) - D(t ) = Q~ 1 BT h {t) + Q-'ETTIIShWPhB - S(t)B}. (4.5) 

lim || /°° e~ M F h (t)dt\\ < \im \\(P h Tl - Tl) f e~ Xt S h (t)Bdt\\ 
h-* o "Jo h—*o Jo 

< UmlKP^n-n^Aji^ii 

h—+ o 

< lim ||Phn||(E(A; Aho)B — R( A; Ao)S]| 

A — +0 

+ ||W^/)ni2(A;A 0 )5|| 

+ |jnJl_H(^(A; A 0 )B - R(\] A^BW 
= 0, 

since Ph converges strongly to the identity and 

lim ||fl(A; A^B - R( A; A 0 )B\\ = 0 

h —*0 

for all A [18, Theorem 4.2]. The sequence will now be shown to be equicontinuous . 

The core C is dense in D(>1) which is dense in X. Let be an orthonormal basis 

for R m and for any 6 > 0 define a bounded operator Gs : #"* — ► X by 

Gf(ei) = Xi where G C and H-Be^ — a:*|| < 6/m. 

Then Gs is a bounded linear operator with range contained in C and with 

\\Gs-B\\<6. 

Now H-A/u.G'iH is bounded by some real Mi since it converges to \\A 0 G S \\. Therefore, 

\\F h (t + 8)-F h (t)\\ < HOT -n|| ||/ 0 5 5 h (f + u)i fco Bdu|| 

< ||n|| ||jj 5 fc (f + u)i^du||+ 

< ||n|| \\[S h {t + 8) - S h (i)](B - G 6 )\\. 
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Let M be the constant bounding §h(t) for small h (Equation 5.2). Then 

|| F h (t + 8)- Ffc(t)|| < \\U\\[8MM 1 + 2 MS] 
and for sufficiently small H and any e > 0 we can choose 8 so that 

\\F h (t + 8)-F h (t)\\<e for all h < ff. 

Thus F h is an equicontinuous family and since F h (t) is also clearly uniformly bounded for 
all t and h , it follows from [16, Lemma 2.11] that 

lim sup || Fh(t) - F(t)|| = 0. 

hr-* 0 0<<<r 

Since II and Q~ 1 B* are bounded operators it follows from Equation (4.5) and Theorem 3.3 
that 

lim sup || D h (t) - D(t)|| = 0. 

0 0<t<T 

Finally, the remainder of the proof of convergence of D h {t) to D(t) follows as for the numer- 
ators (4.3). 

Convergence in the graph topology follows from Theorem 4.1. The conclusions of the 
theorem then follow from Theorem 4.2. 

The result is independent of the technique used for stabilization, and does not require 
any assumptions on the stabilizability or detectability of the approximations. For instance, 
let ( N h ,D h ) be a ^-stable r.c.f. of an approximating plant and choose X,Y G M(A a ) such 

that 

XN h + YD h = /. 

Then if the original plant is close enough to the approximation so that, for some r.c.f of the 
original plant, (JV,D), 

N-N h \\ 1 ( 4i6 ) 

d-DhW* \\XY\\’ 

then, any controller designed using the approximating plant will stabilize the original plant. 
Furthermore, the closed-loop responses will be close in norm, and the system will be internally 
CT-stable. The difficult part, of course, is determining whether (4.6) is satisfied. 

5. Conclusions 

The above results validate the use of Galerkin type approximations m controller design 
for infinite-dimensional systems, regardless of the control technique, and with very weak 
assumptions on the approximation method. The key assumption is that of stabilizability/ 
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detectability of the semigroup control system which supports the work of Jacobson and Nett 
[14] indicating the importance of this property in control of infinite-dimensional systems. 

A high order Galerkin model can be chosen and used to design a controller using lumped 
parameter techniques. If the objective is optimal state feedback, then convergence of the 
Kalman gains can be used as a criterion. Otherwise, robustness of the approximate plant with 
respect to the original infinite- dimensional system can be estimated before any controller is 
implemented by either (1) comparison of coprime factors (if available) or (2) determination 
that the sequence of approximations has converged satisfactorily. Several possible criteria are 
given in [21] and [11]. An open problem is determining when the approximation order is high 
enough for satisfactory controller design. Since exact transfer functions will be unavailable 
for most practical problems, a useful result would be the rate of convergence of the coprime 
factors, compared to that of the numerical scheme. 

A subsequent paper will extend the results in this paper to discrete time approximation 
schemes. 
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